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]^ ' Abstract 

I Continuous-time autoregressive moving average (CARMA) processes have recently been 

' used widely in the modeling of non-uniformly spaced data and as a tool for dealing with 

high-frequency data of the form FnAj'T- = 0, 1,2, . . ., where A is small and positive. Such 
data occur in many fields of application, particularly in finance and the study of turbulence. 
This paper is concerned with the characteristics of the process (K„A)n6Z) when A is small 
and the underlying continuous-time process {Yt)t£M is a specified CARMA process. 
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■ 1 Introduction 

' Throughout this paper we shall be concerned with a CARMA process driven by a second-order 

zero-mean Levy process L with ELi = and ELl = . The process is defined as follows. 

For non-negative integers p and q such that q < p, a. CARMA(p,q) process Y = {Yt)t£R, 
with coefficients oi, . . . , Op, b^, . . . ,bq G R, and driving Levy process L, is defined to be a strictly 
^ , stationary solution of the suitably interpreted formal equation, 

a{D)Yt = b{D)DLt, t G M, (1.1) 

where D denotes differentiation with respect to t, a(-) and 6(-) are the polynomials, 

a{z) := + aizP~^ + ■ ■ ■ + Op and b{z) := 6o -|- biz + •••-!- bp-iz^~^, 

and the coefficients bj satisfy bq = 1 and bj = for q < j < P- The polynomials a(-) and b{-) are 
assumed to have no common zeroes, and it will be assumed that the zeroes of the polynomial a 
all lie in the interior of the left half of the complex plane. 
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Since the derivative DLt does not exist in the usual sense, we interpret (jl.ip as being equiv- 
alent to the observation and state equations 



dX.t = AX.tdt + GpdLt , 



(1.2) 
(1.3) 



where 
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and A = — oi for p = 1. 



A, which we shall denote by Ai, Ap, are 



It is easy to check that the eigenvalues of the matrix 

the same as the zeroes of the autoregressive polynon iial a(-). 

Under the conditions specified it has been shown ( Brockwell and Lindner ( 20091 ) . Lemma 2.3) 
that these equations have the unique strictly stationary solution, 



Y, 



g{t - u)dLu, 



where 



sit) 



zt 



a{z) 



0, 



if t > 0, 



if t < 0. 
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and p is any simple closed curve in the open left half of the complex plane encircling the zeroes 
of a(-). The sum is over the distinct zeroes A of a(-) and ReSz=\{-) denotes the residue at A of 
the function in parentheses. Evaluating these residues, we can write g more explicitly as 

1 



= E 



(m(A) - 1)! 1 



2?r^^)-i [{z - A)"We^*6(z)/a(z) _^ l(o,oo)(t), (1-6) 



where m(A) denotes the multiplicity of th e zero A and D^, denote s diffe rentiation with respect 
to z. The kernel g can also be expressed ( Brockwell and Lindner ( 20091 ) . equations (2.10) and 
(3.7)) as 

<7(t) = bTe^Sl(o,oo)(t)- (1-7) 
From this equation we see at once that g is infinitely differentiable on (0, 00) with /c*'^ derivative, 

gi.k){t) = b^e^U^ep, < t < 00. 

Since bq = 1 and bj = for j > q, the right derivatives 5^^^(0+) satisfy 



^W(0+) = bT^% 



if k < p — q — 1, 

1 if k = p — q — 1, 



fl.8) 
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and in particular g{0+) = 1 ii p — q = 1 and g{0+) = if p — q > 1. 

Gaussian CARMA processes, of whic h the Gaussi an Ornstei n-Uhlenbeck process is an early 
example, were first studied in detail by DoobI (1944) (see also Doob ( 1990)). T he state-space 
formulation, (1.2) and (1.3) (with b"^ = [1 ••• 0]) was used by IjonesI (jl98ll ) to carry out 
inference for time series with irregularly-spaced observations. This form ulation leads natu rally 
to the definition of Levy-driven and non-linear CARMA processes (see iBrockwelll (j200ll ) and 
the references therein). Fractioii ally integrated Levy-driven CARMA processes were studied by 
Brockwell and Marquardtl (j2005l ). 

Levy-driven CARM A processes have been applied successfully to the mo delli ng of stochastic 
volati lity in finance (see lTodorov and TauchenI (j2006l ). lBrockwell et al.l (120061) andlHaug and Czado 



2007)), extending the celebrated Ornstein-Uhlenbeck model of iBarndorff-Nielsen and Shephard 



200 ll ) . The resul ts presented here were motivated by preliminary studies of high-frequency tur- 



bulence data (see iFerrazzand (j20ld )) which appear to be well-fitted by a continuous time moving 
average process sampled at times 0, A, 2 A, . . where A is small and positive. We return to this 
topic in iBrockwell. Ferrazzano and KliiDpelberd (l201lh. The a pplication to turbulence data will 
be investigated in detail in iFerrazzano and Kliippelbergi ()201ll ). 

Our paper is organised as follows. In Section [2] we derive an expression for the spectral density 
of the sampled sequence := (l^A)nGZ- It is known that the filtered process {(j){B)Yj^)n£Z, 
where 4>{B) is the filter defined in (3.1), is a moving average of order at most p — 1. In Section[3l 
we determine the asymptotic behaviour of the spectral density and autocovariance function of 
{(l){B)Y^)n^z as A 4, and the asymptotic moving average coefficients and white noise variance 
in the cases p — q = 1,2 and 3. In general we show that for small enough A the order of the 
moving average {(f){B)Y^)n£Z is p — 1. 



2 The spectral density of := {YnA)nez 

From (jl.Sp we immediately see, since g{t) = for t < 0, that the Fourier transform of g is 



~g{u) := / g{t)e^'^'dt 



b{—iuj) 

27ri J p a{z) z + iui a[—iuj) ' 



h{z) 



dz 



(2.1) 



Since the autocovariance function 7y(-) is the convolution of crg(-) and ag{—-), its Fourier trans- 
form is given by 



7Y(a;) = a g{uj)g{-uj) = a 



b{iu}) 



a{iuj) 



The spectral density of Y is the inverse Fourier transform of 7y . Thus 



Substituting this expression into the relation 

^Y{h) = [ e"^^fY{uj)du;, h G I 



b{iuj) 



a{iLo) 



and changing the variable of integration from u: to z = iu gives, 

.2 r u.\u_.\ . _ n{z)h{-z) 

—-— r-e' • uz = (J y ntitiz=x I 

^ ^ n( 'y\n{ — r 1 ^ — ^ \ 



2m Jp a{z)a{-z) ^ 



a{z)a{—z) 



(2.2) 
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where the sum is again over the distinct zeroes of a(-). 

We can now compute the spectral density of the sampled sequence := (l^A)nGZ- This 
sp ectral density /a will play a key role in the subsequent analysis. We have, from Corollary 4.3.2 
in 



Brockwell and David (Il99ll ). 



2tt 

h=—oo 



and, substituting for 7y from (j2.2p . 

-a^ f b{z)b{-z) sinh(Az) 

/a(w) = ^^/^^7 r — TTT-^ r^dz, -vr < w < vr. (2.3) 

47r^z Jp a[z)a[—z) cosh(A2:j — cos(a;j 

3 The filtered sequence, {(f){B)Yj^)neZ 



If Ai , ■ . . , Ap are the (not necessarily distinct) zeroes of a(-), then we know from lBrockwell and Lindner 
(|2009l ). Lemma 2.1, that if we apply the filter 

cP{B):=fl{l-e^^^B) (3.1) 
i=i 

to the sampled sequence, 

yA 

, we obtain a strictly stat i onary sequence which is {p — l)-correlated 
and is hence, by Lemma 3.2.1 of lBrockwell and DavisI (1199 ll ). a moving average process of order 
p — 1 or less. 

Our goal in this section is to study the asymptotic properties, as A J, 0, of the moving average 
9{B)Zn in the ARMA representation, 

4>{B)Y^ = 9{B)Zn, nGZ, (3.2) 

of the high-frequency sequence V^. Here B denotes the backward shift operator and (^n)nez is 
an uncorrelated sequence of zero-mean random variables with constant variance which we shall 
denote by r^. 

We shall denote by /ma the spectral density of {6{B)Zn)n£Z- Then, observing that the power 
transfer function of the filter (13.11) is 



V;(a;) = | J](l-e^^-^+^'^)|2 = 2Pe-'^i^JJ(cosh(AiA)-cos(t^)), -vr < a; < tt, (3.3) 

j=i i=i 

we have 

/ma(w) = '0(^^)/a(w), -vr < o; < tt, (3.4) 

where 1^(0}) and f^{oj) are given by (j3.3p and (j2.3p respectively. 

In principle the expression (|3.4p determines the second order properties of {9{B)Zn)n& and 
in particular the autocovariances jMAih) for h = 0, . . . ,p — 1. Ideally we would like to use these 
autocovariances to find the coefficients 61, . . . ,9p^i and white noise variance r^, all of which 
are uniquely determined by the autocovariances, if we impose the condition that 6{-) has no 
zeros in the interior of the unit circle. Determination of these quantities is equivalent to finding 
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the corresponding factorization of the spectral density /ma (see ISayed and KailathI (j200ll ) for a 
recent paper on spectral factorization). 

From (j2.3p . (|3.3p and (j3.4p we can calculate the spectral density fMA{^) as —a'^7p(uj)/{2Tr) 
times the sum of the residues in the left half plane of the integrand in (I2.3P , i.e. 



fMAiu.) = -^Hu.) V DT^^y^ ( sinh(A.)b(z)6(-z) _\ 

2T,^y>Z^ i (cosh(AA)-cosH)a(-z)n^^A(^-/^r^''^ ' 



2 = A 



where the sum is over the distinct zeroes A of a(-) and the product in the denominator is over 
the distinct zeroes /i of a(-), which are different from A. The multiplicities of the zeroes A and /i 
are denoted by m(A) and m{^) respectively. When the zeroes Ai, . . . , Ap each have multiplicity 
1, the expression for fMA{^) simplifies to 

(-2re-"^^a^ ^ b{X,)b{-x,) . ^^^^ ^^-rr^ 

JMA{UJ) = > , -77T-r~7 — r^smh(AiA) (cosn^ - cosh(AjA)), -vr < a; < vr. 

Itt ^ a'[Ai)a[-Ai) 

Although in principle the corresponding autocovariances jmaU) could be derived from Jma, 
we derive a more direct explicit expression later as Proposition 13.61 The asymptotic behaviour 
of /ma as a 4, is derived in the following theorem by expanding (j2.3p in powers of A and 
evaluating the corresponding coefficients. Here and in all that follows we shall use the notation, 
a(A) ~ 6(A), to mean that limA4.o a(A)/6(A) = 1. 

Theorem 3.1. The spectral density Jma of {9{B)Zn)n& in the ARMA representation (13. 2p of 
the sampled process has the asymptotic form, as A J, 0, 
2 

fMAico) ~ ^(-l)^'-'^-iA2(f-«)-icp_g„i(a;)2P-i(l - cosu;)^ -vr < c < ^, (3.6) 

where Ck{^) is the coefficient of x^^'^^ in the power series expansion 

. , 00 

Y,c,{u:)x^'^\ (3.7) 



cosh X — cos uj 



k=0 



In particular, co{uj) = j^^, ci(a;) = - gffcos^)^ ' ^2(0;) = ^^'^SolT-otZf^^ 

Proof. The integrand in (12. 3p can be expanded as a power series in A using (13.71). The integral can 
then b e evaluated term by term using the identities, (see Example 3.1.2.3. of 
(|l984h ) 

1 r ^....^^ ^ J , . . {0. 1. 2. . . .,. 

2tti Jp a[z)a[—z) 2 \ a[z)a[—z) J 

from which we obtain, in particular, 

1 /•^2fc+i K^)K-z) /O if 0< A:<p-g-l, 



Mitrinovic and Keckic 



2m Jp a{z)a{-z) \ Ll)LA k = p-q-l. 



Substituting the resulting expansion of the integral (|2.3p and the asymptotic expression 
ip{uj) ~ 2^(1 — cosw)^ into (|3.4p and retaining only the dominant power of A as A — )■ 0, we 
arrive at (13.611. □ 
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Corollary 3.2. The following special cases are of particular interest. 



p-q = l 



p-q 



p-q 



/ma(w) 

fMA{uj) 
fhlAi^) ' 



27r 

271 



2^(1 - cosa;)'5 



2 cos uj\ „, 



cosw 



fj^A^ /II 13costJ cos(2a;)\ 



27r 120 



30 



60 



(3.8) 
(3.9) 
(3.10) 



Proof. These expressions are obtained from (|3.6p using the values of Co(f^), ci{uj) and 02(1^) given 
in the statement of the theorem. □ 

Remark 3.3. (i) The right-hand side of (j3.8p is the spectral density of a g-times differenced 
white noise with variance cr^A. It follows that, ii q = p — 1, then the moving average polynomial 
0{B) in (j3.2p is asymptotically (1 — B)'^ and the white noise variance is asymptotically fi^A 
as A —7- 0. This result is stated with the corresponding results for p — q = 2 and p — g = 3 in the 
following corollary. 



(ii) By Proposition 3.32 of iMarquardt and Stelzeii (j2007l ) a CARMA(p, g)-process has sample 
paths which are {p — q — l)-times differentiable. Consequently to represent processes with non- 
differentiable sample-paths it is necessary to restrict attention to the case p — q = 1. It is widely 
believed that sample-paths with more than two derivatives are t oo smooth to repres e nt th e 
processes observed empirically in finance and turbulence (see e.g. iJacod and Todorovl (j2010l ). 
■Tacod et al.l toid )) so we are not concerned with the cases when p — q > 3. □ 



Corollary 3.4. The moving average process '■= 9{B)Zn in (13. 2p has for A 4, the following 
asymptotic form. 



Xn = il- BfZn, n e 



(a) If p — q = 1, then 

where t'^ := Var{Zn) = a^A. 
(h) Ifp-q = 2, then 

Xn = {l + eB){l-ByZn, neZ, 
where 9 = 2 - y/3 and t"^ := Var{Zn) = a'^A^{2 + V3)/6. 
(c) If p — q = 3, then 

Xn = {l + 9iB + e2B^){l-ByZn, nGZ, 



where 6*2 = 2 (8 + V30) - \/375 + 64^30, Oi = 266*2/(1 + 62) = 13 - \/l35 + 4^30 and 
2 (8 + ^30) + \/375 + 64^/30) A^a'^/120. 



Proof, (a) follows immediately from Theorem 4.4.2 of lBrockwell and DavisI (jl99ll ). 

To establish (b) we observe from (j3.9p that the required moving average is the q times 
differenced MA(1) process with autocovariances at lags zero and one, 7(0) = 2(T^A'^/3 and 
7(1) = fJ^A^/6. Expressing these covariances in terms of 9 and gives the equations, 



{1 + 9^)t^ = 2a^A^/3, 
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from which we obtain a quadratic equation for 6. Choosing the unique solution which makes the 
MA(1) process invertible gives the required result. 

The proof of (c) is analogous. The corresponding argument yields a quartic equation for 62- 
The particular solution given in the statement of (b) is the one which satisfies the condition that 
9{z) is nonzero for all complex z such that |z| < 1. □ 

Although the absence of the moving- average coefficients, 6j, from Corollary 3.4 suggests that 
they cannot be estimated from very closely-spaced observations, the coefficients do appear if 
the expansions are taken to higher order in A. The apparent weak dependence of the sampled 
sequence on the moving-average coefficients as A J, is compensated by the increasing number 
of available observations. 

In principle the autocovariance function ^ma can be calculated, as indicated earlier, from 
the corresponding spectral density Jma given by (13. 4p and (12. ip . Below we derive a more direct 
representation of 7ma and use it to prove Theorem 13.71 which is the time-domain analogue of 
Theorem 13.11 

Define B/^g{t) = g{t - A) for t G M. We show that 0(-Ba)5(-) = for t > pA. 

Lemma 3.5. Let Y be the CARMA(p, g) process (fLi|) and A > 0. Define <i){B) as in ([3T]) . 
Then 

p 

cj){BMt):=ll{l-e^^^BA)g{t)=0, t > pA. (3.11) 

j=i 

Proof. Rewriting the product in (13. IIP as a sum we find (j){B^)g{t) = X]j=o -^jSi^ ~ J^)' which 
has Fourier transform (invoking the shift property and the right hand side of (|2.ip ) 



a{—iijj) 



where the product is taken over the distinct zeroes of a(-) having multiplicity m[\). Using the 
fact that the product of Fourier transforms corresponds to the convolution of functions, we obtain 
from <^ 



Now note that, for every of the distinct zeroes Aj, 

lini ^-^—^ 1^ = A™(^^). 

z^\^ {z - Xj)'"(^j) 

The singularities at z = Xj are removable and, therefore, using Cauchy's residue theorem. The- 
orem 1 of Section 3.1.1, p. 25, and Theorem 2 of Section 2.1.2, p. 7, of iMitrinovic and Keckic 
(jl984l ). the filtered kernel is zero for every i G M. □ 
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Proposition 3.6. Let Y be the CARMA{p,q) process (II. 4p and A > 0. The autocovariance at 
lag n of {(t){B)Y^^)j(,i is, for n = 0,1, ... ,p - 1, 

p—n n+i— 1 1 /.jA 

7MA(n) = ^'Y1 EE / 9{s - hA)gis - {k - n)A)ds, (3.12) 

i=l fc=0 h=0 ^{i-l)A 

Al = {-!)'' e^(^'i+-+^»J, A: = l,...,p. (3.13) 

{n,...,ifc}GC^ 

T/ie SMTO m (|3.13p is taken over the (^) subsets of size k of {1,2, . . . ,p}. 

Proof. We note that ^MAi^) is the same as ^[{(p{Bi\)Y)t{(p{Bi\)Y)t-^-An] and use the same ex- 
pansion as in the proof of Lemma I3.5| i.e. 

<I>{Ba) = flil - b^^^Ba) = E^^i^i, (3.14) 

j=l k=0 

which we apply to Y. Observe that for t G M, setting := t — kA for k = 0,...,p, and 
tp+i ■= -oo, 

BiYt = Bi / g{t - u)dLu = / g{tk - u)dL^ = / gitk - u)dL,,. (3.15) 

J —oo J —CO Jti+i 

Applying the operator (|3.14p to If, using ()3.15p and interchanging the order of summation gives 

A)Y)t=Y. J2^l9(.tk-n)dLu. (3.16) 

From Lemma 13.51 we know that the contribution from the term corresponding to m = p is zero. 
By stationarity, the autocovariance function is independent of t, hence we can choose t = An. 
Then we obtain 

P-l /■A(n-j) i 

{HBA)Y)nA = E / E ^laiin - k)A - u)dLu 

P /■A(n-j+l) i-1 

= E / Y.^l9{{n-k)A-u)dK. 
j=i -^^{"-i) k=o 



For t = 0, we obtain analogously 

P A{i-1) »-l 

iHBA)Y)o = Y, E Kai-^h - u)dLu. 

i=l -^-Ai 

For the autocovariance function we obtain for n = 0, . . . ,p — 1 by using the fact that L has 
orthogonal increments, 

7MA(n) = m(l>iBA)Y)o{(l>iBA)Y)nA] 

p-n i+n-1 i-1 — {«-l)A 

= E E E / . ^((^ - - ^)9i-^h - u)du. 

1=1 A:=0 h=0 J-tA 

Finally, (|3.12p is obtained by changing the variable of integration from ti to s = —u. □ 
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Theorem 3.7. The autocovariance function ^MAin) for n = 1, . . . ,p — 1 has for A ^ the 
asymptotic form 



2a2(p— g)— 1 P""- rt+i~l i-1 / \ / \ 



i=l fc=0 /i=0 

where for iV G No 



C{h,k,i,n;N) : = 



/ (s + i-l-/i)^(s + i-l-A; + n)^ds. 

^0 

Proof. We can rewrite the integral in (|3.12p as 

a[ g{{s + i-l-h)A)g{{s + i-l-k + n)A)ds. (3.18) 
Jo 

Since 5 is infinitely differentiable on (0, 00) and the right derivatives at exist, the integrand has 
one-sided Taylor expansions of all orders M S N, 



[g{{s + i-l- h)A)g{{s + i-l-k + n)A)] 



M 

E dA^ 

1=0 

M I 



A 



-jr + oiA 

A=0+ 



= EE(^P + ^- l- ^)'""(" + i-l-k + nrg('-^^ (0+)5("^) (0+) — + (A^O , 

as A 4- 0. Choose M = 2{p — q — 1). Then by (jl.Sp there is only one term in the double sum 
which does not vanish, namely the term for which m = p — q — 1 = I — m. Setting N := p — q—1 
(so that M = 2N) the sum reduces to 

(^^^ys + ^-l-hfis + ^-l-k + nfJ^A^^ + o{A^^). 

Since (^^) /i2N)] = (iV!)~^, the integral in (fXT8|) is for A | asymptotically equal to 

/ {s + i-l-hf{s + i-l-k + nfds + o{A^^+^), 
Jo 



m 

and, since 



(3.19) 



{ii,...,ih}eCl 

we also have 



hm y e''^^n+-+K)=(P). 



AlAl = {-ir^'(P^I^^+o{l) as A 10. (3.20) 
Combining (fXT9]l and dSJO]), we obtain IKTT} . □ 
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Remark 3.8. (i) For computations the following expansion may be useful (as usual we set 
0° = !) 

C{h,k,i,n;N) := [ {s + i - 1 - h)^ {s + i - 1 - k + n)^ ds 
Jo 



h,l2=0 V 1/ V 2/ Jo 



V ff 1 -, ^--r{^ - 1 - hf-^^{^ -l-k + n^^"'^ 



Furthermore, we observe that C depends on p and q only through p — q. 

(ii) Note that the right hand sides of ()3.5p and ()3.6p are the discrete Fourier transforms of (j3.12p 
and (|3.17p . respectively. Note also the symmetry between (|3.5p and (|3.12p in the dependence on 
A and p — q — 1. □ 

So far we know that the moving average process X„, = 9{B)Zn from (|3.2p is of order not 
greater thanp— 1 but possibly lower. Our next result presents an asymptotic formula for ")ma{p— 
1), which shows clearly that this term is not 0. 

Corollary 3.9. For lag n = p — 1 the autocovariance formula (j3.17p reduces to 

2^2(p-g)-l 

^■«'''-')~'-"' (2(p-,-l)). P'^" 

and ^ma{p — 1) is therefore non-zero for all sufficiently small A > 0. 
Proof. From the expansion p.l7p we find 

Set d := p — q > 1, then 



^^^^P - 1) ~ (l) ^(0' kA,p-l;p-q- 1). (3.22) 



C{0,k,l,p-l;d-l)= [ s'^-^{s-k + p-lf~'ds, 

Jo 

and, from Remark 13.81 this is a polynomial of order d—1. In order to apply known results on the 
difference operator, we define the polynomial f{x) = s'^~^(x + s + p — l)'^~^ds. Then, using 



Eq. (5.40), p. 188, and the last formula on p. 189 in lGraham et al.l (|1994| ). the sum in (|3.22p can 
be written as 

E(-l)'0C(O,A:,l,p-l;d-l) 

= - ^)l-=o - (-1)'' Q ^(0' k,l,p-l;d-l) 

= + (-l)P+i r s'^~\s - if-^ds = (-1)^+"^ s'^-\l - sf-^ds, (3.23) 
Jo Jo 
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p-q 


1 


2 


3 


4 


1ma{p - 1) 


A(-1)P-V2 


6-1a3(-1)p-V2 


120-iA5(-l)P-3a2 


5040-iA^(-l)P-V^ 



Table 1: Values of 7a/a(p — 1) for p — q = 1, . . . , 4. 



where we have used the fact that d — 1 = p — q — 1 < p. To obtain Eq. (|3.2ip it suffices to note 
that = {-lfP-1 = (-1)'' and that the integral in lK23\i is a beta function. Hence 



T{2d) (2d-l)! 

□ 

Remark 3.10. If Y is the CARMA(p, q) process (1.4) then, from Theorem 13.11 the spectral 
density of (1 — B)p~'^Y^ is asymptotically, as A J, 0, 

2 

^A(-2A2)P-9-icp_g_i(a;)(l - cosw)P-«, vr < a; < vr. 

If p — 5 = 1, 2 or 3 this reduces to the corresponding spectral densities in Corollary 3.2, each 
divided by 2*^(1 — cosw)"^. The corresponding moving average representations are as in Corollary 
3.4 without the factors (1 - B)i . 

In particular, for the CAR(l) process, (1 — B)Y^ has a spectral density which is asymptot- 
ically cr^A/(27r) so that, in the Gaussian case, the increments of Y^ for small A approximate 
those of Brownian motion with variance a'^t. □ 



In this pape r we have considered only second-order properties of y^. It is possible (see 



Brockwelll ( 200ll ). Theorem 2.2) to express the joint characteristic functions, E exp(i X^^^ ^k^h^ 



k 



for m G N, in terms of the coefficients Oj and hj and the function ^(•), where ^{6) for ^ € M 
is the exponent in the characteristic function Ee*^^^ = e^^^^ of Li. In particular the marginal 
characteristic function is given by Eexp(iey;f) 

These expressions are awkward to use in practice, however Brockwell. Davis and Yan^ ( 201ll ) 
have found that least squares estimation (which depends only on second-order properties) for 
closely and uniformly spaced observations of a CARMA(2,1) process on a fixed interval [0, T] 
gives good results. They find in simulations that for large T the empirically-determined sample 
covariance matrix of the estimators of ai , 02 and 69 is close to the matrix calculated from the 
asymptotic (as T — t- 00) covariance matrix of the maximum likelihood estimators based on 
continuous observation on [0, T] of the corresponding Gaussian CARMA process. 



4 Conclusions 

When a CARMA(p, q) process Y is sampled at times nA for n € Z, it is well-known that the 
sampled process Y^ satisfies discrete-time ARMA equations of the form (j3.2p . The determination 
of the moving average coefficients and white noise variance for given grid size A, however, is a 
non-trivial procedure. In this paper we have focussed on high frequency sampling of Y . We have 
determined the relevant second order quantities, the spectral density fMA of the moving average 
on the right-hand side of (|3.2p and its asymptotic representation as A J, 0. This includes the 
moving average coefficients as well as the variance of the innovations. We also derived an explicit 
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expression for the autocovariance function ^ma and its asymptotic representation as A J, 0. This 
shows, in particular, that the moving average is of order p — 1 for A sufficiently small. 
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